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METHOD FOR CHARACTERISING AN OPTIC AL FIBRE LINK 

10 Field of the invention 

[0001] The present invention is related to a method 

for characterising an optical fibre link by analysis of a 
Rayleigh backscatter signal, from which the spatial 
distribution of the beat length, the coupling length and 

15 the PMD (Polarization Mode Dispersion) along the link can 
be determined. 



State of the art 

[0002] The use of dispersion-shifted and dispersion 

20 compensating fibres has minimised the effect of chromatic 
dispersion on the bandwidth of an optical link. 
Polarisation mode dispersion (PMD) has therefore become the 
most serious limiting factor in high-speed optical 
communication systems. There exist several measurement 
25 techniques of PMD, but they only allow for measuring the 
global value and do not give information about the 
distribution of PMD along the fibre length. Measuring the 
spatial distribution of PMD is, however, important. Such a 
measurement would indeed allow to locate the bad trunks in 
30 an optical link, causing high PMD value. This can be 
essential in the frame of a network maintenance and 
upgrade . 

[0003] The PMD of a fibre depends on two parameters: 

the beat length (Lb) and the coupling length (L c ) . The beat 
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length depends on the birefringence of the fibre and the 
coupling length is related to the phenomenon of 
polarisation mode coupling. Measuring the spatial 
distribution of PMD therefore consists of measuring the 
5 distribution of both the beat and coupling lengths. 

[0004] Most of the techniques proposed for 

measurement of the distributed PMD are based on 
polarisation-optical time domain ref lectometry (POTDR) . The 
concept of POTDR was introduced by A. Rogers in the 80 's 
10 (see A. Rogers, * Polarization-optical time domain 
ref lectometry: a technique for the measurement of field 
distributions', Applied Optics, vol. 20, pp. 1060-1074, 
1981) . It basically consists of measuring the polarisation 
properties of the Rayleigh backscatter signal when an 
15 optical pulse propagates down the fibre. Some POTDR-based 
techniques have since been developed for PMD 
characterisation. Most of them allow a measurement of the 
beat length distribution along an optical link and, 
therefore, do not allow the complete determination of the 
20 PMD. However, Sunnerud presented a POTDR-based set-up 
enabling the measurement of the accumulation of PMD along 
the fibre (see H. Sunnerud et al., "Polarization -Mode 
Dispersion measurements along installed optical fibers 
using gated backscattered light and a polarimeter' , J. of 
25 Lightwave Technol., vol. 18, pp. 897-904, 2000). This 
method requires POTDR measurements at several wavelengths. 
A technique based on the measurement of the degree of 
polarisation by using a POTDR set-up has also been 
presented in "Distributed detection of high-PMD sections on 
30 installed fibers using a polarization-OTDR' , M. Leblanc, 
Proc. OFMC'01, pp- 155-162, 2001 , and is related to the 
patent application US 2003/0174312. This technique does not 
allow quantification of the PMD, but enables to obtain 
distributed information about the level of PMD. Moreover 
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this technique requires the complete measurement of the 
state of polarisation of the backscattered signal. More 
recently, Galtarossa described another technique for the 
PMD measurement for which the coupling length is determined 
5 from the correlation of the birefringence vector (see 
^Measurement of birefringence correlation length in long, 
single -mode fibers ' , Gal tarossa et al. , Opt. Let t . , vol . 
26, pp. 962-964, 2001) , which also requires the measurement 
of the polarisation states of the backscattered signal. 

10 Moreover, this technique requires a quite small duration of 
the pulses launched in the optical fibre: 5-10 ns. 
[0005] In x Polarization mode dispersion mapping in 

optical fibers with a Polarization-OTDR' , AT. Wuilpart, G. 
Ravet, P. Megrret and M. Blondel, Photon. Technol. Lett., 

15 vol. 14, pp. 1716-1718, 2002, the authors describe an 
analysis of the POTDR signal which allows a mapping of the 
PMD along an optical fibre link by quantifying the PMD on 
each fibre of the link. The main advantage of that 
technique is that it does not require the complete 

2 0 measurement of the backscattered polarisation state 
evolution and only uses a linear polariser at the fibre 
input, which is quite simple to implement. The 
determination of the beat length, the coupling length and 
finally the PMD is based on the analysis of the statistical 

25 properties of the extrema present in a POTDR trace. 

[0006] ' In said paper by Wuilpart et al. the spatial 

distribution of the linear birefringence S is characterised 
by a Rayleigh statistical distribution along the fibre 
length and it is assumed that the rate of change of the 

30 birefringence angle q is described by a Gaussian 
distribution with zero mean and standard deviation a 
depending on the polarisation mode coupling strength. The 
circular birefringence is neglected. The corresponding 



measurement set-up is shown in Fig.l. The OTDR pulses (10 
ns) externally modulate a 1553 nm DFB laser via a pulse 
generator and an electro-optic modulator (EOM) . After 
amplification, the pulses are launched into the fibre 
through an acousto-optic modulator (AOM) , which suppresses 
the amplified spontaneous emission noise of the erbium 
doped fibre amplifier between two successive pulses. A 
linear polariser is placed at the fibre input and a 
polarisation controller (CP) is used to obtain the maximum 
power after the polariser. 

[0007] Fig. 2 shows three simulated POTDR traces 

corresponding to the POTDR arrangement and obtained for 
different values of the beat length and an identical <7 
value. One can easily observe that the number of maxima 

strongly depends on Lb- 

[0008] The calculation of the beat length is based 

on the measurement of the number of maxima of the POTDR 
trace. Fig. 3, which has been obtained by simulation, 
presents the number n of maxima per unit of length in 
function of the beat length for different values of a. This 
figure shows that n varies not only with L B but also with 
the strength of the polarisation mode coupling (a) . In 
order to estimate lis from n, the curve corresponding to a = 
10 (the dashed curve in Fig. 3) is used. This choice indeed 
minimises the maximum possible error if one assumes that a 
is smaller than 15 degrees/m. Consequently, this method 
allows only to give an approximate value of the beat length 
L B . A relationship between the coupling length L c and cr can 
be derived and one can obtain: 

0.879 

L C ~ 2~~ 

cr 

[0009] A statistical approach is then used for 

determining the coupling length from the POTDR trace. Fig. 
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4 shows different POTDR traces that were obtained by 
simulation for a beat length of 25 m and for different 
values of a: 2, 4 and 8 degrees/m, respectively. One can 
observe that the lower envelope of the curves varies in a 
5 smoother way when a decreases and therefore when L c 
increases. Defining APi as the absolute value of the 
difference in ordinate between two successive minima of the 
POTDR trace P one can write : 

A^=|P(z,)-P(z #+I )| 

10 where Zi and z i+1 correspond to the distances of the i th and 
the i+1 minimum, respectively. 

[0010] A new parameter § can then be calculated as 

£ =< AP. >. 

where < . >i denotes the mean value for i varying from 1 to 
15 the number of minima of the POTDR trace. £ is a measure of 
the speed variation of the lower envelope of the 
backseat tered signal along the fibre length. Fig. 5 
(obtained by simulation) shows the evolution of £ according 
to a for different values of the beat length L B - The values 

20 are here expressed in dB because £ has been calculated 
directly from the POTDR trace. 

[0011] The procedure for the determination of L c is 

the following. From the POTDR trace, the number of maxima 

# 

and each AP± are measured. An approximate value of the beat 
25 length is then deduced by the method previously described 
and £ is calculated by averaging APi. Finally a can be 
deduced from £ by using the curve of Fig. 5 corresponding 
to the correct beat length and then Lc is calculated by 
means of the relation Ii C =0.879/a 2 stated above. Because L B 
30 is an approximate value, the calculated coupling length L c 
will also be an approximation. 



[0012] * From the measurement of the beat length and 

the coupling length, it is easy to determine an approximate 
value of the PMD value by 



PMD 




\2 f „ 2 ^ 



5 where c is the light velocity in a vacuum, X the 
measurement wavelength and L the fibre length. It should be 
emphasised that the final PMD value is approximate: this is 
due to the approximation done on the beat length 
measurement . This process can be repeated for each 
10 contribution in the POTDR signal corresponding to the 
different fibres of the optical link: it is therefore 
possible to map the PMD along this link. 

[0013] The method as described above has some 

important drawbacks and limitations. It only yields 

15 approximate values of Lb, L c and the PMD. Further, the 
presented POTDR traces were assumed ideal. It supposes 
ideal characteristics of the various components of the 
experimental set-up. In practice, the set-up is 
characterised by some imperfections, which cause a 

20 distortion of the ideal POTDR signal. Therefore the method 
as such cannot be directly applied. The POTDR signal is 
also affected by noise: detector noise and residual 
coherence noise, which adds a series of minima and maxima 
to the POTDR trace. The final value of the measured PMD can 

25 therefore be erroneous. 



Aims of the invention 

[0014] The present invention aims to provide a 

method for analysing a backscatter signal for determining 
30 in an accurate way the spatial distribution of the beat 
length, the coupling length and the polarisation mode 
dispersion. 
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Summary of the invention 

[0015] The invention relates to a method for 

characterising an optical fibre link by its beat length, 
5 coupling length and polarisation mode dispersion 
distribution. It comprises the steps of 

- sending a pulsed signal along the optical fibre link and 
measuring the backscattered signal (preferably a POTDR 
signal) , after passing through a polariser, 
X0 - deriving the length of the optical fibre, the average 
power difference between two successive minima of the 
backscattered signal and the number of maxima per unit 
length, 

in an iterative way determining a beat length interval 
15 and an interval foar the polarisation mode coupling 

parameter, until the length of the intervals is below a 
predetermined value, yielding a value for the beat 
length and the coupling length, 

calculating the polarisation mode dispersion. 
20 [0016] Preferably the POTDR signal is an ideal POTDR 

signal . 

[0017] Advantageously the POTDR signal is the 

convolution of an ideal POTDR signal and a signal depending 
on the pulse shape. In an alternative embodiment the POTDR 
25 signal further is convoluted with a signal taking into 
account the time jitter introduced by the measurement set- 
up. 

[0018] In a specific embodiment a smoothing 

algorithm is applied to the POTDR signal. 
30 [0019] In another object the invention relates to a 

method for characterising an optical link consisting of a 
concatenation of several fibres, wherein the method as 
previously described is applied to each fibre. 
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[0020] The invention also relates to the use of a 

method as described above to locate the position of 
polarisation mode dispersion sources within an optical 
fibre link. 

[0021] In an advantageous embodiment the method as 

described is used in telecommunication networks. 
Alternatively, it can be applied in fibre sensing 
applications . 

Short description of the drawings 

[0022] Fig. 1 represents the set-up for measuring a 

POTDR trace. 

[0023] Fig. 2 represents three simulated POTDR 

traces for a given a and for different Ii B values. 

[0024] Fig. 3 represents the number of maxima per 

unit length versus L B for different values of cr . 

[002 5] Fig. 4 represents simulated traces for a 

given beat length and different values of <j. 

[0026] Fig. 5 represents the evolution of £ as a 

function of a for different values of Ii B . 

[0027] Fig. 6 represents the determination of the 

beat length interval . 

[0028] Fig. 7 represents the determination of the 

interval for a. 

[0029] Fig. 8 represents the dis torsion of the POTDR 

signal due to the imperfections of the set-up (pulse shape 
+ time jitter) . 

[0030] Fig. 9 represents the evolution of n with the 

beat length when the imperfections of the set-up are taken 
into account . 

[0031] Fig. 10 represents the evolution of £ with a 

when the imperfections of the set-up are taken into 
account . 



[0032] Fig. 11 represents the simulated and measured 

probability density function of A. 

[0033] Fig. 12 represents the window algorithm that 

can be used for calculated the local PMD. 



[0034] It is first shown how the polarisation mode 

dispersion from a POTDR signal can be determined in an 
accurate way CLb is not approximately determined anymore) 

10 by analysing the statistical properties of its extrema. 
Suppose first that a falls within the range [0,15] 
degrees/m, which means that the coupling length is assumed 
to be larger than 12-8 m. This <y interval has been chosen 
such that it corresponds to the values measured on real 

15 fibres. 

[0035] From a POTDR signal, the fibre length L, £ 

and the number of maxima per unit of length, n*, are easily 
measurable. After measuring n^, it is possible to determine 

a first interval \^B l0W > ^B up ] for the beat length of the 
20 fibre under test so that: 



r^ 1 >l i i 



B ^ and B up are determined by the intersections of the 

curves n(L B ,a) corresponding to a = 0 and a = 15 degrees/m 
with the straight line n = rim as shown on Fig. s. 
25 [0036] When the first beat length interval is 

determined, the curves ^(a,L B ) can be used to define an 

interval [^/ou/^ip] for the a value of the fibre under test 
so that : 



<J e |cr /( 
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cr} ow and °up are determined by the intersections of the 



1 r 7I 

curves £(ct,Lb) for L B= L B and L B =L B with the straight 

tow 



line J; = £ m# where £ m is the £ value of the fibre under test 
measured from the POTDR signal- This principle is 
illustrated in Fig. 7. 

[0037] After the determination of <?} 0 w and &l p , the 

curves n(IiB/tf) for = and & = &up can be calculated- A 

new beat length interval is then deduced as previously 
explained by Fig. 6, but using the curves obtained for 

cr = <j} ow and °' = 0 'lp instead of a = 0 and a = 15 degrees /m. 
[0038] This new beat length interval denoted 

^ ow >^J is included in [^!ff low ,jt ^] and gives a smaller 
range for the beat length of the fibre under test. Fig. 7 

[ 2 2 1 

is then used again to calculate a new a interval \Fhw> cr up\ 
by using the curves £(cy,L B ) obtained for 

L B = L 2 R and L B = L \ 

This new interval gives a smaller range for the <j value of 
the fibre under test. This process is repeated until small 
enough intervals are obtained for the beat length L B and 

for <y. 

[0039] When the beat length Lb and a are determined, 

the PMD is finally calculated by means of the relationship: 



PMD = J— 
V 3tt 



i 



1.545 



n 



K L B a Cj 



tr 2 L 



0.439 



— 1 + e 0 439 



[0040] If the optical link to be characterised 

consists of a concatenation of several fibres, the method 
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is applied to each fibre of the link and their beat and 
coupling lengths are therefore determined independently of 
each other. In this way, the PMD mapping is possible along 
the link. 

5 [0041] The method for analysing the backscatter 

signal as described offers several advantages : 

- it requires only the access to one end of the optical 
fibre link to be characterised; 

- it is based on a POTDR arrangement, but does not require 
10 the complete measurement of the complete state of 

polarisation of the backscattered signal. Only a simple 
polariser is required ; 

it involves a POTDR measurement at only one wavelength; 

- it gives accurate (not approximate) values for the beat 
15 length, the coupling length and the PMD. 

[0042] In practice, the final measured signal is 

also affected by noise. This noise is due to the detector 
noise and to the presence of a residual coherence noise. As 
the determination of the PMD is based on the detection of 

■ 

20 extrema, it is obvious that this noise will affect the 
results . 

[0043] As already mentioned before, the measured 

signal is also affected by the imperfections of the 
experimental set-up like the finite pulse shape and the 

25 time jitter of the electronic equipment. 

These phenomena (noise and imperfections of the measurement 
set-up) lead to a distortion of the POTDR signal. Therefore 
it cannot be used directly to deduce the PMD map. 
The width of the optical pulses launched into the fibre 

30 under test is finite and can be non-rectangular. 
Consequently, the backscattered power measured at the POTDR 
corresponds to an average of the ideal signal over a small 
spatial window depending on the pulse shape. The resulting 



12 



signal measured at the POTDR is in fact the convolution 
product between the ideal signal and the shape of the 
pulse : 

PpulseW = nz)®Pshape{z) 

where P pu ise(z) is the power measured by the POTDR according 
to a certain pulse shape. P(z) is the ideal POTDR signal 
and P S hape(z) is the power distribution along the pulse. 
Hence the Fourier transform of P pu ise{z) becomes; 

P pulsed = n^PshapeW 

where P pu ise(k), P(k) and P 3 ha P e(k) are the Fourier transforms 
of PpuiseU), P(z) and P shaP e(z), respectively. 

[0044] Moreover the electronic equipment is not 

perfect: it is subject to a time jitter, i.e. the time 
delay between two successive pulses is not constant. This 
phenomenon strongly affects the averaging process of the 
POTDR. If we suppose that the number of pulses involved in 
the averaging process of the OTDR is large enough and that 
the time jitter is a random variable characterised by a 
Gaussian distribution of zero mean and standard deviation 
<y t j, the signal detected at the POTDR can be written as: 



P m W = PpulseW®e « 

where P ra (z) is the power measured by the POTDR taking into 
account the pulse shape and the time jitter. Hence the 
Fourier transform of P m (z) becomes: 

P m (k) = P pulse (k)a tj e 2 

The pulse shape and the time jitter have for effect to 
distort the measured backscattered signal . An example is 
shown on Fig. 8, which has been obtained for parameters of 
the experimental set-up. The signals (a) and (b) correspond 
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to the ideal and affected signals, respectively. One can 
clearly observe that the signal is strongly distorted. 
[0045] Because of these imperfections, the method 

described above cannot be used directly: the curves n(L B ,a) 
and ^ (or, Lb) are indeed based on ideal POTDR signals 
obtained by simulations. Therefore the effects of the pulse 
shape and the time jitter should be included in the 
calculation of n(Le,a) and E,(o,L B ) . In the case of the 
parameters of the experimental set-up and a pulse duration 
of 10 ns, these new curves are shown on Fig. 9 and 10. In 
order to deduce these curves, ideal POTDR traces have been 
simulated and the effects of the imperfections have been 
taken into account by using the above expressions for 
Ppuise(z) and P a (z) . Consequently, when the imperfections of 
the experimental set-up are taken into account, the PMD map 
of an optical link can still be determined.. The same 
reasoning as the one explained previously has to be applied 

i 

on the POTDR signal by using the new set of curves n(L B ,cr) 
and £(ct,L b ) described in fig. 9 and 10. As the pulse shape 
is now taken into account, it is possible to launch a 
larger pulse in the fibre (50 ns, for example), which will 
increase the system dynamics. The maximum measurable length 
will therefore be greater. 

[0046] As already indicated, the final measured 

signal is also affected by noise. This noise is due to the 
detector noise and the presence of a residual coherence 
noise. Its main effect is to add a series of minima and 
maxima to the POTDR trace. The values of n and £ cannot 
thus be correctly measured and the final calculated PMD 
will be erroneous. In order to solve this problem, a 
smoothing algorithm is applied to the backseat tered POTDR 
signal P m (z). This algorithm consists of generating the 
signal P B (z) such that: 
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P s (z) is therefore the mean value of P m (z) on a distance 
interval L 9 around z. The number of maxima and the 
distances z k corresponding to the locations of the minima 
of the POTDR trace are then calculated from the smoothed 
signal in order to ignore the influence of the small peaks 
generated by the noise. Afterwards, when the positions z k 
are known, % is calculated from P m (z) • When noise is 
present, the number of extrema of P s (z) depends of the 
smoothing parameter L a . So, the main problem is to choose 
the correct smoothing parameter to be applied to the POTDR 
signal . 

[0047] The location of the minima of the smoothed 

signal can be used to estimate the probability density 
function of A, denoting the difference in z between two 
successive minima. This estimation can then be compared to 
the density probability function of A of a POTDR trace 
obtained by means of simulations taking into account the 
parameters Lb and a calculated from the smoothed signal. If 
the probability density functions of A obtained from P 3 (z) 
and from the simulation do not match, it means that the 
noise contribution is significant. The smoothing parameter 
L 3 will therefore be chosen such that the corresponding 
probability density function obtained from measurements can 
be correctly fitted by the probability density function 
derived from simulations: L s is therefore chosen such that 
the noise contribution is reduced. Fig. 11 illustrates this 
concept. Figures (a), (b> and (c) show the probability 
density functions obtained by measurement (circles) and by 
simulation (crosses) for a given optical fibre and for an 
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value of 3, 4 and 5 m, respectively. One can clearly 
observe in this case L s = 3 m gives the best fitting. 
[0048] The invention also provides a technique which 

will allow to locate more accurately the position of the 
PMD source within the fibre. The basic idea consists in 
applying the described method on a distance window (1 km, 
for example), which slides along the fibre length. By 
comparing the obtained PMD's on several successive windows 
as shown on Fig. 12, it is possible to locate more 
accurately high PMD sources within the optical fibre. 
[0049] POTDR can also be used in fibre sensing. The 

polarisation properties indeed depend on several external 
effects like strain, temperature, electric and magnetic 
fields. Strain, twist, temperature, electric and magnetic 
fields affect the birefringence within the fibre and 
consequently, the mean beat length. The method of the 
invention, which measures the spatial distribution of the 
beat length, therefore allows to measure the spatial 
distribution of these external effects after calibrating 
the sensing fibre. 
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CLAIMS 

1. Method for characterising an optical 
link by its beat length, coupling length and 
polarisation mode dispersion distribution, comprising the 
steps of 

- sending a pulsed signal along said optical fibre link 
and measuring the backseat tered signal, after passing 
through a polariser, 

- deriving the length of said optical fibre, the average 
power difference between two successive minima of said 
backscattered signal and the number of maxima per unit 
length, 

- in an iterative way determining a beat length interval 
and an interval for the polarisation mode coupling 
parameter, until the length of said intervals is below a 
predetermined value, yielding a value for the beat 
length and the coupling length, 

- calculating the polarisation mode dispersion. 

2. The method as in claim 1, wherein said 
backscattered signal is a POTDR signal. 

3. The method as in claim 2, wherein said 
POTDR signal is an ideal POTDR signal. 

4. The method as in claim 2, wherein said 
POTDR signal is the convolution of an ideal POTDR signal 
and a signal depending on the pulse shape. 

5. The method as in claim 4, wherein said 
POTDR signal further is convoluted with a signal taking 
into account the effect of time jitter. 

6. The method as in any of claims 2 to 5, 
wherein a smoothing algorithm is applied to said POTDR 
signal . 



7. Method for characterising an optical link 
consisting of a concatenation of several fibres, by 
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applying the method as in any of the previous claims to 
each fibre. 

8. Use of method as in any of the previous 
claims to locate the position of polarisation mode 
dispersion sources within an optical fibre link. 

9. Use of method as in any of the claims 1 
to 8 in telecommunication networks. 

10. Use of method as in any of the claims 1 
to 8 in fibre sensing applications. 
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The present invention is related to a method 
for characterising an optical link by its beat length, 
coupling length and polarisation mode dispersion, 
comprising the steps of : 

- sending a pulsed signal along said optical fibre link 
and measuring the backscattered signal (being a POTDR 
signal) , after passing through a polariser, 

- deriving the length of said optical fibre, the average 
power difference between two successive minima of said 
backscattered signal and the number of maxima per unit 

5 length, 

- in an iterative way determining a beat length interval 
and an interval for the polarisation mode coupling 
parameter, until the length of said intervals is below a 
predetermined value, yielding a value for the beat 

0 length and the coupling length, 

- calculating the polarisation mode dispersion. 
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